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BETHE ANSATZ AND THE SPECTRAL THEORY OF AFFINE 
LIE ALGEBRA-VALUED CONNECTIONS. 

THE SIMPLY-LACED CASE 

DAVIDE MASOERO, ANDREA RAIMONDO, DANIELE VALERI 


Abstract. We study the ODE/IM correspondence for ODE associated to 
g-valued connections, for a simply-laced Lie algebra g. We prove that subdom¬ 
inant solutions to the ODE defined in different fundamental representations 
satisfy a set of quadratic equations called ’L-system. This allows us to show 
that the generalized spectral determinants satisfy the Bethe Ansatz equations. 
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0. Introduction 

The ODE/IM correspondence, developed since the seminal papers [13, 3], con¬ 
cerns the relations between the generalized spectral problems of linear Ordinary 
Differential Equations (ODE) and two dimensional Integrable Quantum Eield The¬ 
ories (QET), or continuous limit of Quantum Integrable Systems. Accordingly, the 
acronym IM usually refers to Integrable Models (or Integrals of Motion). More re¬ 
cently, it was observed that the correspondence can also be interpreted as a relation 
between Classical and Quantum Integrable Systems [27], or as an instance of the 
Langlands duality [17]. 

The original ODE/IM correspondence states that the spectral determinant of 
certain Schrodinger operators coincides with the vacuum eigenvalue Q{E) of the Q- 
operators of the quantum KdV equation ]4], when E is interpreted as a parameter 
of the theory. This correspondence has been proved ]13, 3] by showing that the 
spectral determinant and the eigenvalues of the Q operator solve the same functional 
equation, namely the Bethe Ansatz, and they have the same analytic properties. 

The correspondence has soon been generalized to higher eigenvalues of the Quan¬ 
tum KdV equation - obtained modifying the potential of the Schrodinger operator 
]5] - as well as to Conformal Field Theory with extended >V„-symmetries, by consid¬ 
ering ODE naturally associated with the simple Lie algebra A„ (scalar linear ODEs 
of order n-l-1) ]10, 34]. Within this approach, the original construction corresponds 
to the Lie algebra Ai, and the associated ODE is the Schrodinger equation. 
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After these results, it became clear that a more general picture including all 
simple Lie algebras - as well as massive (i.e. non conformal) Integrable QFT - 
was missing. While the latter issue was addressed in [20, 27] (and more recently in 
[11, 1, 30[), the ODE/IM correspondence for simple Lie algebras different from A„ 
has not yet been fully understood. 

The first work proposing an extension of the ODE/IM correspondence to clas¬ 
sical Lie algebras other than is [9[, where the authors introduce a g-related 
scalar ODE as well as the important concept of d'-system. The latter is a set of 
quadratic relations that solutions of the g-related ODE are expected to satisfy in 
different representations. It was then recognized in [17] that the g-related scalar 
differential operators studied in [5, 9[ can be regarded as ajfine opers associated to 
the Langlands dual of the affine Lie algebra g. This discovery gave the theory a 
solid algebraic footing - based on Kac-Moody algebras and Drinfeld-Sokolov Lax 
operators - which was further investigated in [32[ and that we hereby follow. 

In the present paper we establish the ODE/IM correspondence for a simply laced 
Lie algebra g, in which case the affine Lie algebra g coincides with its Langlands dual 
[17[. More concretely, we prove that for any simple Lie algebra g of ADE type there 
exists a family of spectral problems, which are encoded by entire functions 
i = 1,... ,n = rankg, satisfying the ADE-Bethe Ansatz equations [31, 35, 2[: 

(n^E*) 

= -1 > ( 0 - 1 ) 

for every E* G C such that Q^^'>{E*) = 0. In equation (0.1), C = is 

the Cartan matrix of the algebra g, while the phase D and the numbers /3j are the 
free-parameters of the equation. More precisely, following [17[ we introduce the 
g-valued connection 

C{x,E) = dx + - + e+p{x,E)eo , (0.2) 

X 

where £ € t) is a generic element of the Cartan subalgebra 1) of g, eo, ei,..., Cn are the 
positive Chevalley generators of the affine Kac-Moody algebra g and e = 

In addition, the potential^ is p{x, E) = x^^'^ — E, where M > 0 and h'^ is the dual 
Coxeter number of g. 

After [32[, for every fundamental representation ^ of g we consider the 
differential equation 

C{x,E)'^{x,E) =0. (0.3) 

The above equation has two singular points, namely a regular singularity in x = 0 
and an irregular singularity in x = oo, and a natural connection problem arises 
when one tries to understand the behavior of the solutions globally. We show that 
for every representation there exists a unique solution ^'(*)(x, E) to equation 
(0.3) which is subdominant for x —>■ -l-oo (by subdominant we mean the solution 
that goes to zero the most rapidly). Then, we define the generalized spectral 
determinant Q^'^'^{E;i) as the coefficient of the most singular (yet algebraic) term 
of the expansion of 4'^*)(x, E) around x = 0. Finally, we prove that for generic £ € t) 
the generalized spectral determinants Q^^\E;£) are entire functions and satisfy the 
Bethe Ansatz equation (0.1), also known as Q-system. 

The paper is organized as follows. In Section 1 we review some basic facts 
about the theory of finite dimensional Lie algebras, affine Kac-Moody algebras and 

^We stick to a potential of this form for simplicity. However all proofs work with minor 
modifications for a more general potential discussed in [5, 17]. 

^Actually yb) is an evaluation representation of the i-th fundamental representation of g. 
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their finite dimensional representations. Furthermore, we introduce the relevant 
representations that we will consider throughout the paper. 

Section 2 is devoted to the asymptotic analysis of equation (0.3). The main result 
is provided by Theorem 2.4, from which the existence of the fundamental (subdom¬ 
inant) solution 4'’^*)(x, £1) follows. The asymptotic properties of the solution turn 
out to depend on the spectrum of A = cq -|- e G g, as observed in [32] . 

The main result of Section 3 is Theorem 3.6, which establishes the following 
quadratic relation among the functions d'(®^(x,£i), known as 'h-system, that was 
conjectured in [9]: 

{x,E) A ■^f{x,E)) = (x, £;)®^« . (0.4) 

^2 2 ' 

Here, 4'^\(x,£i) is the twisting of the solution 'h(®^(x,£i) defined in (2.8), while 

rui is the morphism of g-modules defined by (1.11), and B = denotes 

the incidence matrix of g. In order to prove Theorem 3.6, we first establish in 
Proposition 3.4 some important properties of the eigenvalues of A on the relevant 
representations previously introduced. More precisely, we show that in each repre¬ 
sentation there exists a maximal positive eigenvalue (see Definition 2.3), 
and that the following remarkable identity holds: 

n 

(^e~^ + i = 1,..., n = rankg . (0.5) 

Equation (0.5) shows that the vector ..., A*^"^) is the Perron-Frobenius eigen¬ 
vector of the incidence matrix B, and this implies that the A^*^’s are (proportional 
to) the masses of the Affine Toda Field Theory with the same underlying Lie algebra 
0 [18]. 

In Section 4 we derive the Bethe Ansatz equations (0.1). To this aim, we study 
the local behavior of equations (0.3) close to the Fuchsian singularity x = 0, and we 
define the generalized spectral determinants Q^'^\E;£) and Q^''\E;i) using some 
properties of the Weyl group of g. In addition, using the ^'-system (0.4) we prove 
Theorem 4.3, which gives a set of quadratic relations among the generalized spec¬ 
tral determinants (equation (4.6)). Evaluating these relations at the zeros of the 
functions Q^''\E]i), we obtain the Bethe Ansatz equations (0.1). 

Finally, in Section 5, we briefly study an integral representation of the sub¬ 
dominant solution of equation (0.3) with a linear potential p(x, E) = x, and we 
compute the spectral determinants Q^''\E]tj. Some explicit computations of the 
eigenvectors of A are provided in Appendix A. 

In this work we do not study the asymptotic behavior of the {E) ’s, for E ^ 0, 

that was conjectured in [9], neither we address the problem of solving the Bethe 
Ansatz equations (0.1) - for a given asymptotic behavior of the functions Q^''\EYs 
- via the Destri-de Vega integral equation [8]. These problems are of a rather 
different mathematical nature with respect to the ones treated in this work, and 
they will be considered in a separate paper. 


Note added in press. We proved the ODE/IM correspondence for non simply- 
laced Lie algebras in [29]. 
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1. Affine Kac-Moody algebras and their finite dimensional 

REPRESENTATIONS 

In this section we review some basic facts about simple Lie algebras, affine Kac- 
Moody algebras and their finite dimensional representations which will be used 
throughout the paper. The discussion is restricted to simple Lie algebras of ADE 
type, but most of the results hold also in the non-simply laced case; we refer to 
[19, 24] for further details. At the end of the section we introduce a class of g-valued 
connections and the associated differential equations. 

1.1. Simple lie algebras and fundamental representations. Let g be a simple 
finite dimensional Lie algebra of ADE type, and let n be its rank. The Dynkin 
diagrams associated to these algebras are given in Table 1. Let C = 
be the Gartan matrix of g, and let us denote by R = 21„ — C the corresponding 
incidence matrix. Since g is simply-laced, it follows that Bij = 1 if the nodes i,j 
of the Dynkin diagram of g are linked by an edge, and Bij = 0 otherwise. 
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Table 1: Dynkin diagrams of simple Lie algebras of ADE type. 


4 


Let {fi,hi,ei \ i G I = n}} C g be the set of Chevalley generators of g. 

They satisfy the following relations {i,j G /): 

hj] — 0 , [hit ; [hit fj\ — fj ; [^ii fj] — ^ijhi , (1.1) 

together with the Serre’s identities. Recall that 

f) = Chi C g 
iei 

is a Cartan subalgebra of g with corresponding Cartan decomposition 



As usual, the Qa are the eigenspaces of the adjoint representation and ii C f)* is the 
set of roots. For every i G I, let ai G [}* be defined by ai{hj) = Cij, for every j G I. 
Then Qa, = Cccn- Hence, A = {ai \ i G 1} C R. The elements ai,...,a„ are 
called the simple roots of g. They can be associated to each vertex of the Dyiikin 
diagram as in Table 1. We denote by 

P = {A G r I Kh^) G Z, Vi G /} C 1)* 

the set of weights of g and by 

P+ = {A G P I Xih,) > 0, Vi G /} C P 

the subset of dominant weights of g. We recall that we can define a partial ordering 
on P as follows: for A, /x G P, we say that A<^ifA — ^isa sum of positive roots. 
For every A G (}* there exists a unique irreducible representation P(A) of g such 
that there is u G P(A) \ {0} satisfying 

hiV = \{hi)v , CiV = 0 , for every i G I . 

L(A) is called the highest weight representation of weight A, and v is the highest 
weight vector of the representation. We have that P(A) is finite dimensional if and 
only if A G P’*', and conversely, any irreducible finite dimensional representation of 
g is of the form P(A) for some A G P"*". For A G P'^, the representation P(A) can be 
decomposed in the direct sum of its (finite dimensional) weight spaces P(A)^, with 
/X G P, and we have that P(A)^ 0 if and only if ^ < A. We denote by P\ the set 

of weights appearing in the weight space decomposition of P(A); the multiplicity 
of /X in the representation P(A) is defined as the dimension of the weight space L(A)^. 

Recall that the fundamental weights of g are those elements uji G P"*", i G It satis¬ 
fying 

u!i{hj) = Sij t for every j G I. (1-3) 

The corresponding highest weight representations P(axi), i G It are known as funda¬ 
mental representations of g, and for every i G I we denote by Vi G L{u}i) the highest 
weight vector of the representation L{uji). Since the simple roots and the funda¬ 
mental weights of g are related via the Cartan matrix (which is non-degenerate) by 
the relation 

iOi — ^ ^ i G I t 

then we can associate to the x—th vertex of the Dynkin diagram of g the corre¬ 
sponding fundamental representation L(uji) of g. 


5 


1.2. The representations Lijii). Let us consider the dominant weights 

7^2 - ^ ^ ^ij 5 i ^ I 

j<^I 

where B = is the incidence matrix defined in Section 1.1, and let L{r]i) 

be the corresponding irreducible finite dimensional representations. In addition, we 
consider the tensor product representations 

i€l. 

We now show that we can find a copy of the irreducible representation L{r]i) inside 
the representations /\ L{uJi) and This will be used in Section 3 

to construct the so-called ■^-system. 

Lemma 1.1. The following facts hold in the representation /\^ L{uji). 

(a) For every j G I, we have that 

CjifiVi A Ui) = 0 . 

(b) For every h € (}, we have that 

h{f^Vi A Vi) = rii{h){f^Vi A Vi). 

(c) Any other weight of the representation /\^ L{uJi) is smaller than rji. 

Proof. Since Vi is a highest weight vector for L(uJi) we have that CjVi = 0 for every 
j G L Moreover, using the relations [ej,fi] = Sijhi together with equation (1.3), 
we have 

V-jfiVi — fiV-jVi F SijhiVi — dijVi , 

and therefore ej{fiVi A Vi) = SijVi A Vi = 0, proving part (a). By linearity, it 
suffices to show part (b) for h = hi, i G L Recall that, by equation (1.1), we have 
[hi,fj] = -Cijfj, for every i,j G L Then, 

[^jfi\) Cji)fiVi . 

Hence, hj{fiVi A u*) = {2Sij - Cji){fiVi A Vi) = r]i{hj){fiVi A Vi). Finally, let w ^ 
0 Ji,r]i — uji be a weight appearing in LitOi). It follows from representation theory of 
simple Lie algebras that 

ui < r]i - uJi = uJi - CjiUjj < uJi. 

Therefore, by part (b), the maximal weight of /\^ L{uji) is Wi + rji — uJi = rji, thus 
proving part (c). □ 

As a consequence of the above Lemma it follows that fiViAvi G /\^ L{uJi) is a highest 
weight vector of weight rji. Since g is simple, the subrepresentation of /\ L[(jJi) 
generated by the highest weight vector fiVi A Vi is irreducible. Therefore, it is 
isomorphic to L(rji). By the complete reducibility of L{iOi) we can decompose 
it as follows: 

2 

l\L{ujf) = L{rj^) ®U , 

where U is the direct sum of all the irreducible representations different from L{rji). 
By an abuse of notation we denote with the same symbol the representation L(rji) 
and its copy in /\^ L{uji). For every i G I, let us denote by Wi = The 

following analogue of Lemma 1.1 in the case of the representation (Slje/ 
holds true. 

Lemma 1.2. The following facts hold in the representation . 
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(a) For every j G I, we have that 


CjWi = 0. 

(b) For every /i G (}, we have that 

hwi = r]i{h)wi. 

(c) Any other weight of the representation is smaller than rji. 

Proof. Same as the proof of Lemma 1.1. □ 

As for the previous discussion, by Lemma 1.2 it follows that Wi € L{ujj)^^*^ 

is a highest weight vector which generates an irreducible subrepresentation isomor¬ 
phic to L{rji). Lemmas 1.1 and 1.2, together with the Schur Lemma, imply that 
there exists a unique morphism of representations 

2 

mi : f\L{u}i) ^ , (1.4) 

3^1 

such that Ker mi = U and mi{fiVi A vf) = Wi. 

1.3. Affine Kac-Moody algebras and finite dimensional representations. 

Let hP be the dual Coxeter number of g. Let us denote by k the Killing form of g 
and let us fix the following non-degenerate symmetric invariant bilinear form on g 
(a,5 G g): 

ia\b) = ^K{a\b). 

Let g 0 be the space of Laurent polynomials with coefficients in g. For 

a{t) = S we let 

Rest=o a{t)dt = a_i. 

The affine Kac-Moody algebra g is the vector space g = g 0 C[t, t~^] © Cc endowed 
with the following Lie algebra structure (a, & G g, f{t), g{t) G C[t, 

[a®f{t),b®g{t)\ = [a,b](^f{f)g{t) + {a\b)Rest=o{f'{f)g{t)dt)c, 

[c.a = o. 

The set of Chevalley generators of g is obtained by adding to the Chevalley genera¬ 
tors of g some new generators /o, ho, eo (for the construction, see for example [24]). 
The generator cq, which plays an important role in the paper, can be constructed 
as follows. There exists a root —9 G i?, known as the lowest root of g, such that 
—9 — ai ^ R, for every i € I. Then 

eo = a © t, for some a G g_e . (1-6) 

We now consider an important class of representations of g. Let K be a finite 
dimensional representation of g. For C G C* we define a representation of g, which 
we denote by K(C), as follows: as a vector space we take K(C) = V, and the action 
of g is defined by 

{a'Si f{t))v = f{(){av), cv = 0, ior a € g , f{t) € C[t,t~^] ,v € V . 

The representation K(C) is called an evaluation representation of g. If V and W 
are representations of g then any morphism of representations f : V ^ W can 
obviously be extended to a morphism of representations / : K(C) —>■ kk(C)j which 
we denote by the same letter by an abuse of notation. Similarly, when referring 
to weights or weight vectors of the evaluation representation K(C), we mean the 
weights and weight vectors of the representation V with respect to the action of g. 
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For a representation F of g and k € C, we denote Vk = the evaluation 

representation of g corresponding to ^ . Clearly, if A: G Z then 14 = Vb = 

1^(1), and if fc G i + Z then 14 = 14 = 1^(~1)- Then, for each vertex i G / of the 
Dynkin diagram of g, we associate the following finite dimensional representation 
of g: 

=L(w,)p(i) , (1.7) 

2 

where Liuji) is the fundamental representation of g with highest weight and the 
function p : I ^ 'Ljl'L is defined inductively as follows: 

p(l) = 0 , p{i) = p{j) + 1, for j < i such that By > 0. (1.8) 

By an abuse of language we call the representations (1.7), for i G /, the finite 
dimensional fundamental representations of g. These representations will play a 
fundamental role in the present paper. 

Example 1.3. In type An, = L(a;i)o is the evaluation representation at C = 1 
of the standard representation Lipji) = Moreover, 

i 

for fc= l,...,n. 

' ' 2 

Example 1.4. In type £)„, = L{u!i)o is the evaluation representation at C = 1 

of the standard representation L{uji) = C^". Moreover, 

i 

yW = A yli) , for A: = 1,..., n - 2 . 

' ' 2 

Finally, = L{(jJn-i)^ and = L{u}n)^ are the evaluation representations 

at ( = (—1)" of the so-called half-spin representations. 

Example 1.5. In type Eq, = L(wi)o and are the evaluation 

at C = 1 of the two 27-dimensional representations (they are dual to each other). 
Moreover, 

1/(2) ^ yy yW ^ y(3) ^ yy yC) ^ ^ y{4) ^ ^ y{5) 

Finally, E(®( = L{ujq)i is the evaluation representation at ^ = — 1 of the adjoint 
representation. 

As a final remark, note that by equation (1.7) we have 

(1-9) 

V / P(0+1 ^ 

2 

and if we introduce the evaluation representation 

M(*( =(g) (l/(4)®-®«, VzGl, (1.10) 

then by equation (1.7), and from the fact that By ^ 0 implies p{i) = p{j) -\- 1, it 
follows that the morphism rrii defined by equation (1.4) extends to a morphism of 
evaluation representations 

2 

rui : A ^ . (1.11) 

' ' 2 

The construction of the '0—system will be provided by means of the above extended 
morphism. Due to Lemmas 1.1 and 1.2, the evaluation representation 

1F« =B(77,,)p(0+i 4g J, 

2 
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( 1 . 12 ) 




is a subrepresentation of both and M^'^\ and the copies of Vbb) jn these 

2 

representations are identified by means of the morphism rrii. 

1.4. g-valued connections and differential eqnations. Let {fi,hi,ei \ i = 
0,..., n} C g be the set of Chevalley generators of g, and let us denote e = 

Fix an element £ G i). Following [17] (see also [32]), we consider the g-valued 
connection 

£{x, E) = dx + — + e + p{x, E)eo , (1-13) 

X 

where p{x, E) = x^^ — E, with M > 0 and E G C. Since e is a principal nilpotent 
element, it follows from the Jacobson-Morozov Theorem together with equation 
(1.6) that there exists an element h G £) such that 

[h,e]=e, [/i,eo] =-(h'^ - l)eo, (1.14) 

see for instance [7[. Under the isomorphism between f) and 1)* provided by the 
Killing form, the element h corresponds to the Weyl vector (half sum of positive 
roots). Let k G C and introduce the quantities 

27ri 2-KiM 

uj = , U = eT^=a;'‘". 

Moreover, let be the automorphism of g-valued connections which fixes dx, g 
and c and sends t -G Then from equations (1.5), (1.13) and (1.14) we get 

Mk{uj’^'^'^'"C{x,E)) =uj'^C{uj'^x,^tE). (1.15) 

We denote Ck{x,E) = Aik {jC-{x, E)), for every k G C. Note in particular that 
Co{x,E) = C{x,E). Equation (1.15) implies that for a given hnite dimensional 
representation V of g, the action of Ck{x,E) on the evaluation representation Vq 
of g is the same as the action of C{x, E) on 14, for every k G C. In addition, let C 
be the universal cover of C*, suppose ip{x, E) : C —>■ Vb is a family - depending on 
the parameter E - oi solutions of the (system of) ODE 

C{x, E)ip{x, E) = 0 , (1-16) 

and for k G C introduce the function 

<Pk{x,E)=uj-^’^p{uj'^x,n^E). (1.17) 

Since on any evaluation representation of g we have 

gada^^ga^g-a, a^b GQ , (1.18) 

then by equations (1.15) and (1.18), we have that 

Ck{x,E)ipk{x,E) =Q, (1-19) 

namely pk{x^ E) : C —>■ Vb is a solution of (1.16) for the representation I 4 . 


2. Asymptotic Analysis 


Let V be an evaluation representation of g. Let us denote by C = C \ M<o, the 
complex plane minus the semi-axis of real non-positive numbers, and let us consider 
a solution : C —>■ U of 


C[x,E)'^{x) = 'I''(a;) 


- -f e+p{x, E)eo 4'(a:) = 0 . 

X 


( 2 . 1 ) 


Equation (2.1) is a (system of) linear ODEs, with a Fuchsian singularity at a; = 0 
and an irregular singularity at a; = 00 . Our goal is to prove the existence of 
solutions of equation (2.1) with a prescribed asymptotic behavior at infinity in a 
Stokes sector of the complex plane. In the present form, however, equation (2.1) 
falls outside the reach of classical results such as the Levinson theorem [26] or the 
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complex Wentzel-Kramers-Brillouin (WKB) method [16], since the most singular 
term p(x, E)eo is nilpotent. In order to study the asymptotic behavior of solutions 
to (2.1) we use a slight modification of a gauge transform originally introduced in 
[32], and we then prove the existence of the desired solution for the new ODE by 
means of a Volterra integral equation. 

We now consider the required transformation. First, we note that the function 
p{x,E)t^ has the asymptotic expansion 

p{x,E)~i^ = q{x,E) + 0{x~^~^), 

where 

^ = M(h'^(l + s)-l)-l>0, (2.2) 

and 

S 

q{x,E) = x^ (2.3) 

For every j = 1,..., s, the function Cj{E) is a monomial of degree j in E. 

Lemma 2.1. Let C{x, E) be the Q-valued connection defined in (1.13) and let h G i) 
satisfy relations (1.14). Then, we have the following gauge transformation 

q{x, E) = dx-\- q{x,E)A + - —, (2.4) 

X 

where A = cq + e, the function q(x, E) is given by (2.3) and S by (2.2). 

Proof. It follows by a straightforward computation using relations (1.14). □ 


Remark 2.2. The transformation (2.4) differs from the one used in [32] and given 
by 


„M ad h 


C{x,E)=da^+x‘^K + 


£-Mh 


CqE 


The latter transformation fails to be the correct if M{h'^ — 1) < 1- In fact, in this 
case the term cqE/x^^^ goes to zero slowly and alters the asymptotic behavior. 


We are then left to analyze the equation 

'^'{x) + ^q{x, E)A+ -—Eh _|_ A(a;)^ 4'(a;) = 0 , (2.5) 

where A{x) = 0(x~^~^^ ) is a certain matrix-valued function. From the general 
theory we expect the asymptotic behavior to depend on the primitive of q{x,E), 
which is known as the action. In the generic case ^ Z_|_, the action is defined 
as 

S{x,E)=f q{y,E)dy, xGC, (2.6) 

Jo 

where we chose the branch of q{x,E) satisfying q ^ \x\^ for x real. In the case 
G Z+ then formula ( 2 . 6 ) becomes 

S{x,E) = '^f Cj{E)y^^^-^''^'idy + Cs{E)\ogx, s = ' 

3=0 “'o 

Existence and uniqueness of solutions to equation (2.1) (or equivalently, of (2.5)) 
depend on the properties of the element A = e-|-eo in the representation considered. 
We therefore introduce the following 

Definition 2.3. Let A be an endomorphism of a vector space V. We say that 
a eigenvalue A of A is maximal if it is real, its algebraic multiplicity is one, and 
A > Re p, for every eigenvalue ^ of A. 
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We are now in the position to state the main result of this section. 


Theorem 2.4. Let V be an evaluation representation of g, and let the matrix 
representing A G 'q in V have a maximal eigenvalue A. Let if & V be the corre¬ 
sponding unique (up to a constant) eigenvector. Then, there exists a unique solution 
'^{x,E) : C —> 1^ to equation (2.1) with the following asymptotic behavior: 

'i>{x,E) = e~^^^^’^\{x,E)~'^(^ijjo{l)) as a: —>■+oo . 


Moreover, the same asymptotic behavior holds in the sector |argx| < 2 (m+i) > 
is, for any 5 > Q it satisfies 

'i'{x,E) = e~^^^^'^\{x,E)~'^(^ijj-\-o{l)) , in the sector | argx| < :r 7 TT—^ — (5. 


2(M+ 1) 


The function 4'(x, E) is an entire function of E. 


(2.7) 


Remark 2.5. We will prove in Section 3 that the hypothesis of the theorem are 

satisfied in the ADE case for the representations V^'^\ ^7^*^ and W^''\ 

2 

which were introduced in the previous section. 


Remark 2.6. In the case M{h'^ — 1) > 1, the asymptotic expansion of the solution 
in the above theorem reduces to 


4'(x)=e ^~^TTTq[x,E) ''(?/> +o(l)) , in the sector | argx| 


< 


2(M+ 1) 


Before proving Theorem 2.4, we apply it to prove the existence of linearly inde¬ 
pendent solutions to equation (2.1). First, note that for any fc € M such that 
|fc| < lL_LM±li^ the function 

^k{x,E)=uj-^^'^{uj'^x,Vl'^E), xeR+ (2.8) 

defines, by analytic continuation, a solution : C —>■ 14 of equation (2.1) in the 
representation 14. Theorem 2.4 implies the following result. 

Corollary 2.7. For any k € M. such that \k\ < IlJM+11^ t/jg positive real 
semi-axis the function has the asymptotic behaviour 

^k{x,E) = e-^^''^^^'^\{x,E)-'^-t-^'^{if + o{l)), x>0, (2.9) 

2-Ki 

where if is defined as in Theorem 2.) and 7 = e . Moreover, let I G Z-|_ be 
such that 1 < ? < ^ . Then, the functions 4* ,..., dt are linearly independent 

solutions to equation (2.1) in the representation Vi+i. 


Proof. A simple computation shows that S{oj^x,Ll^E) = 7^S'(x,A), and similarly 
q{u}^x,Lf E) = uj^q{x,E). Therefore by Theorem 2.4, 

4'(a;'=x, Ft^^E) = uj-^^e-'^'‘^^^^’^^q{x, E)-^{if + o(l)), x » 0 . 


X » 0, 


Hence, on the positive real semi axis we have 

^fe(x, E) = A)-S-"AV' + o(l)), 

where we have used the identity = 7 . Note that 7 ^, k G ■ ■ • i 

are pairwise distinct. Therefore, the vectors ifk are linearly independent. This 
concludes the proof. □ 


Remark 2.8. A slight variation of the proof of Theorem 2.4 presented below allows 
one to prove that the asymptotic behavior of dt holds on a larger sector of the 
complex plane, consisting of three adjacent Stokes sectors. This refined result is 
necessary to consider the lateral connection problem, where one considers the linear 
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relations among solutions that are subdominant in different Stokes sectors. Such a 
problem leads naturally to another instance of the ODE/IM correspondence, namely 
to spectral determinants related to the transfer matrix of integrable systems, see 
e.g. [28], 

2.1. Proof of Theorem 2.4. To complete our analysis, we need to further simplify 
equation (2.5) with a new gauge transformation. Let us denote by C(x,E) = 
q{x,E)'^^C{x,E), where C{x,E) and h are as in Lemma 2.1. 

Lemma 2.9. Let i = and h = '^i^j aihi, with it, at G C, and con¬ 
sider the element N = ~ Mai)fi € 0 . Then we have the following gauge 

transformation: 

e“(^) ^C{x, E)=d + qix, E)A + Oix -^-^), 

where a{x) = {xq{x,E)) ^ . 

Proof. It follows by a straightforward computation using the explicit form of >C(a;, E) 
given in equation (2.4), together with the definitions of N and a{x) and the com¬ 
mutation relations (see ( 1 . 1 )) 

[fi, Go] = 0 and [fi, ej] = S^jhi , for every i,j G I. 

□ 


By Lemma 2.9 and equation (1.18), it follows that the ODE (2.5) is transformed 
into 

$'(x) -b (q{x, E)A + A{x)^ $(a;) = 0 , (2.10) 

where A{x) = 0(x“^“^), and 4'(a:) = G(x)4'(x), with G{x) = Since 

G{x) = 1-1- o{x~^), then we look for a solution to equation (2.10) with the same 
asymptotic behavior (2.7) as the solution 'I'(a;) to equation (2.5). Finally, we define 
$(a;) = Then, it is easy to show that equation (2.10) takes the form 

4>'(x) -b (qix, E)A + A{x)^ $(x) = 0, (2.11) 


where A = A — Al. Note that Aip = 0. The problem is now reduced to prove the 
existence of a solution to equation ( 2 . 11 ) satisfying, for every e > 0 , the asymptotic 
condition 


d'(x) ='ip o(l), 


in the sector | argx| < 


TT 

2 (M-b 1 ) 


— e. 


( 2 . 12 ) 


To proceed with the proof we need to introduce some elements of WKB analysis, 
for which we refer to [16, Chap. 3]. We fix k > 0 (possibly depending on E) big 
enough such that the Stokes sector 


E = {x e C I ReS'(x, E) > Re5 '(k, E)} 


has the following two properties: S' is a bi-holomorphic map from E to the right 
half-plane {zSC|Re 2 ;> 0 }, and E coincides asymptotically with the sector 
I argxj < 2 {m+i) • other words, any ray of argument | argxj < 2 {m+i) eventually 
lies inside E, while any ray of argument [argxj > 2 (m-\-i) eventually lies in the 
complement of E. For any e > 0, we define = S“^{| arg 2 ;| < f — e}, and we 
introduce the space B of holomorphic bounded functions U : —>• R for which the 
limit limaj^oo U (x) = U (oo) is well-defined. The space B is complete with respect 
to the norm ||C/||oo = sup^^g^,^ \U{x)\. 
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We construct the solution of equation (2.11) with the desired asymptotic behavior 
(2.12) as the solution of the following integral equation in the Banach space B: 

<^{x) = K[^]ix) + , (2.13) 

where the Volterra integral operator K is defined as 

K[<i>]ix) = r A{y)^{y)dy . 

J X 

The integral above is computed along any admissible path c, and we say that a 
(piecewise differentiable) parametric curve c is admissible if |c(s)|ds = 0(|c(t)|) 
and Re S{c{t), E) is a non decreasing function. For example, 5'“^(f+a;) is an admis¬ 
sible path [16]. Note that a solution $(x) of equation (2.11), with the asymptotic 
behavior (2.12) belongs to B. Moreover, it satisfies equation (2.10) if and only 
if it satisfies the Volterra integral equation (2.13). Indeed, provided that RT is a 
continuous operator, we have that lima;_>oo Rl[d>](a;) = 0 and 

R'[$]'(a;) = —^(a;)$(a:) — q{x^E)KK[^]{x) 

= —y4(a;)<i)(a:) — q{x, E) A ^d)(a:) — 

= — A(a;)<i)(a:) — q{x, E) A$(a;), 

where the last equality follows from the fact that A'0(®) = 0. 


In order to prove that the solution of the integral equation (2.13) exists and it is 
unique we start showing that AT is a continuous operator on B and its spectral 
radius is zero, i.e. lim„^oo ||Ar”||" = 0. By hypothesis, all eigenvalues of A have 
non-positive real part. Since A is semi-simple, i.e. diagonalizable, we deduce that 
there exists a /3 > 0 such that for any k >0 and any (p € V, 

|e«Vl </3|¥^|. (2.14) 

Since |A(a;)| = 0{x~^~^), we have that |A| G L^{c,\c(t)\dt) for every admissible 
path c , and the function 

pOO 

Pc{x)=^ |A(c(s))||c(s)|ds (2.15) 

J X 

is well-defined and satisfies lim^^-j-oo Pc(a^) = 0. Moreover, the inequality (2.14) 
implies that 



E^i^^^^^^-S(V’E))A{y)U{y)dy 


<Pc{x)\\U\U 


for every U £ B. By (2.14) and (2.15) it follows [16] that K[U] is path-independent 
and thus well defined. Furthermore, we have 


\K[U]{x)\ < \\U\\^p{x) 


and ||A:|| = sup ||AT]!/]jjoo < P, 
|(7|oo = l 


where 

p{x) = inf pc{x ), p = sup p{x). 

c admissible x^Dg 

Note that p(x) is a bounded function such that lima;_>oo p{x) = 0. Hence, we 
conclude that A' is a well-defined bounded operator on B. By definition, AT”]!/] 
can be written in terms of K{x, y) as 


poo poo 

K^[U]{x) = ... K{x,yi)...K{yn-i,yn)U{yn)dyi...dyr, 

J X J Vn — ^ 


Vn-l 
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where K(x, y) is the matrix-valued function 

and from this it follows that K'^[U]{x) can be estimated as 

\K-[U]{x)\ = ||C/||oo r. . r |I(2/l)| . . . |I(yn)Myi . . .dj/n < ||C/||oo^ ■ 

Jx Jyn-1 

Thus 

||K"||<^, (2.16) 

nl 

and the series 

$ = ^ K"[v>] 

neZ-j- 

converges in B. Clearly $ satisfies the integral equation (2.13), since 

OO 

n—1 

Moreover, $ is the unique solution of equation (2.13), for if $ is another solution 
then 

[$ — $] = $ — $, for every n G Z+ . 

By equation (2.16), it follows that $ = l>. 

It remains to prove that $ is an entire function of E. This follows, by a stan¬ 
dard perturbation theory argument, from the fact that A and S are holomorphic 
functions of E in D^. Theorem 2.4 is proved. 

3. The ^'-system 

In this section, for a simple Lie algebra 0 of ADE type, we prove an algebraic 
identity known as 'I'-system, which was conjectured first in [9]. Recall from Section 
1.3 that we can write 
2 

/\ rP ^ IrW © C/, ^IRW©C/, (3.1) 

where the representation = Litji) p(i)+i, i G I, was defined in equation (1.12), 
^ 2 

and U and U are direct sums of all the irreducible representations different from 
IR(*). Due to (1.11), for every i G I we have a morphism of representations 

rrii : A —>■ , Kerm^ = U. 

' ' 2 

The 'k-system is a set of quadratic relations, realized by means of the morphisms 
rrii, among the subdominant solutions to the linear ODE (2.1) defined on the rep¬ 
resentations /\Vi^ and 
2 

In order to prove the main result of this section, we need to study the maximal 

eigenvalue (and the corresponding eigenvector) of A in the representations /\^ Ri*\ 

2 

and IR(®) . To this aim, we recall some further facts from the theory of simple 
Lie algebras and simple Lie groups. First, we need the following 

Lemma 3.1. Let h G i) satisfy relations (1.14), and let us set 7 = e"?^. Then, for 
every k G C-, the following formula holds: 

7'=’"d/*A = 7'=A4_fc(A). 
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(3.2) 




In particular, if 'll) is an eigenvector of A in an evaluation representation V, then 


ifk = eigenvector of A in the representation Vk, and ive have 

Aif = Alp if and onl'y if A'lpk = 'y^Xipk ■ (3-3) 

Proof. It follows directly from the commutation relations (1-14). □ 

Now consider the element A = A|t=i G g, which is well known [25] to be a regular 
semisimple element of g. Therefore, its centralizer g^ = {a G g | [A, a] = 0} is a 
Cartan subalgebra of g, which we denote t)'. Due to Lemma 3.1, it follows that 

7 ’^'^^A = 7 A, (3.4) 

where h G h is the element introduced in Section 1 satisfying relations (1.14). 


Equation (3.4) says that A is an eigenvector of and this in turn implies that 

fl' is stable under the action of Since is a Coxeter-Killing automorphism 

of g [25], we can regard g Coxeter-Killing automorphism of t)'. From these 

considerations it follows - as proved in [ 6 ] - that we can choose a set of Chevalley 
generators of g, say {/',/i',e' | * G /} C g, such that t]' = 0ig/that the 
identity 

21 „+ 7 ad/*_^ 7 -ad?>^(^i) 2 ^ (3.5) 

holds as operators on [)'. Here B denotes, as usual, the incidence matrix of g. 
Following [18], relation (3.5) can be used to express - in the basis {/i',* G /} - 
the eigenvectors of in terms of the eigenvectors of B. Now it is well known 
that the incidence matrix H is a non-negative irreducible matrix [24]. Therefore, 
it has a maximal (in the sense of the Perron-Frobenius theory) eigenvalue, which 
is 72 -|--y “2 j and a corresponding eigenvector, say {xi,X 2 , ■ ■ ■ ,Xn) G R>o- If 
xi = 1, we then have the following relation for every i £ I: 

'^BijXj = + ji)xi, Xi > 0. (3.6) 

jei 

Following [18], and essentially using (3.5), we finally obtain 

^ = (3.7) 

iei 

where Xi,i G I, are defined in (3.6) and the function p{i), i £ I, is given in (1.8). 

Remark 3.2. Note that equation (3.7) corresponds to a precise normalization of A. 
In fact, it is always possible to hnd an automorphism of g, under which f) is stable, in 
such a way that A = cj^i^i 7 ^^ Xih'i for any c G C* and any choice of the relative 

I p(i) 

phase 7 ^“ 2 “. The choice of the relative phase is however immaterial. Indeed, under 

p(0 

the transformation 7 “ 2 “ —)• q 2 “ the right hand side of (3.7) represents a different 
element in the Cartan subalgebra [}' which is however conjugated to A. Hence, their 
spectra coincide. 

The decomposition (3.7) allows one to compute the spectrum of A using the 
Perron-Frobenius eigenvector of the incidence matrix B. We will also need the 
following: 

Lemma 3.3. Let aj, j £ I, be the positive roots of q with respect to the Cartan 
subalgebra . Then A satisfies 

A) = ' (1 — 7 “ 2 p 0 ')-i- 1 ) ^ foralli,j£l, (3.8) 

where the Xi are defined in (3.6). This implies that 

Reaj( 7 "^A) > 0 , if p{i)=p{j), 

Reaj( 7 “^A) = 0, if pii) 7 ^ p{j) ■ 
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(3.9) 


Proof. Since Bij 7 ^ 0 if and only p{i) ^p{j)^ equation (3.6) implies that 

E p(j) _ , _i 1 , 

7 2 = (7 2+72)7 2 Xi . 


jei 


Equation (3.9) and the decomposition (3.7) imply the thesis. 


□ 


Using the above lemma, together with equation (3.7), we compute the maximal 
eigenvalue of A in the representations V^'^\ /\ ]+(*) and 

Proposition 3.4. Let i £ I. If g is a simple Lie algebra of ADE type, then the 
following facts hold. 

(a) For the fundamental representation V^’’^ there exists a unique maximal eigen¬ 
value of A. The eigenvalues € I, satisfy the following identity: 

= ( 7-5 +^^)aW. (3.10) 

3^1 

We denote by '0^*^ S the unique (up to a constant factor) eigenvector of A 
corresponding to the eigenvalue A*-*^ . 

(b) For the representation M^'^\ we have that is a maximal eigenvalue 

of A. Moreover, 

3^1 

is the corresponding eigenvector. 

(c) For the representation /\^ we have that ( 7 “^ + 7 ^)A*-®^ is a maximal eigen- 

2 

value of A. Moreover, 

2 

7 - 5 VW a75'‘0« g a y*’ 

' ' 2 

is the corresponding eigenvector. 

(d) The elements 7 “ 5 ^ 0 b) G V'P and 0^^ G belong to the sub¬ 

representation . Hence, for the representation there exists a maximal 
eigenvalue of A, given by 

pii) = (7”^ +7^)A('\ 

3^1 

and the relation 

TOi( 7 “ 2 '®^ 0 ) 72 '®^!®)^ = 

holds for some c G C*. 

Proof. Recall that the weights appearing in the representation L{uJi) are of the 
following type: 

UJ2: Oij CV , 

where in the latter case the index j is such that 7 ^ 0 and cu is an integral non¬ 
negative linear combination of the positive roots. Clearly for any h € tj' such that 
Reaj(h) > 0, Vj € I we have 

RewAh) > (wi — ai){h) > {uii — ai — Uj — U3){h). 

Therefore for any such h the eigenvalues with maximal real part is uiiih), followed 
by (wi — ai){h), (wi — ai — aj){h), and so on. 

We use this simple argument to study the maximal eigenvalues of A in the 
representation By definition of A and A, and using equation (3.2), the action 
of A on coincides with the action of 7^(“1+®®'^'®)A on L{uji). By Lemma 3.3, 


16 



Q;j(7 ^2'A) > 0, for every j € I, and 0:^(7 ^2'A) = Xi{l — 7^ It follows 

that in the representation V^'^\ A has the maximal eigenvalue 

p(i) — 

Wi(7 5 -A) = Xi. 

Therefore, by definition of the Xi’s, the satisfy equation (3.6). 

Moreover we have 

(3.11) 

and 

= C7-V^*^ , (3.12) 

for some c G C*. This completes the proof of part (a). Part (b) follows directly from 
part (a). To prove part (c), we follow the same lines we used to prove part (a). In 

fact, the action of A on coincides with the action of on L{uji), 

2 

where p*{i) = 1 — p{i). By equation (3.8) we have Reai(7^^“^“^+'^'^^^A) = 0, 

while Reaj (7^^^“^+‘^'^^^A) > 0 if p{i) ^ p{j)- It follows that the two eigenvalues 
with maximal real part correspond to to the weights oji and Wj — a^. By equation 
(3.8), these eigenvalues are Xi^^ and xci~'^^ and the corresponding eigenvectors 
are 7^2 Therefore /\ has a unique maximal eigenvalue ( 72+7 ^)xi, 

with eigenvector 

VW A 75 A /'u'), (3.13) 

for some c ^ 0. Here we have used equations (3.11) and (3.12). 

Let us prove (d). By equation (3.11) and Lemma 1.2, it follows that G 
where IT^®) is seen as a subrepresentation of . Therefore the maximal eigenvalue 
of IT^®) and coincide. Finally by equations (3.11), (3.13) and Lemma 1.1, we 
have that A e and m,( 7 - 5 '‘?/>(b a c G C*, 

thus proving part (d) and concluding the proof. □ 

Remark 3.5. We note that we can always choose a normalization of the eigenvectors 
'0(®) G z G /, such that 

771 ,( 7 -^ A 7 ^ = V'® , (3.14) 

for every z G /. From now on, we will assume that Proposition 3.4(d) holds with 
the normalization constant c = 1 . 

By Proposition 3.4(a), for any fundamental representation V^'‘\ i £ I of g, there 
exists a maximal eigenvalue with eigenvector Therefore, by Theorem 2.4 
there exists a unique subdominant solution '^'^'^\x,E) : C —>■ of the differential 
equation ( 2 . 1 ) in the representation with asymptotic behavior 

4'W(a;,F;) = E)-'^ + o(l)) , in the sector |arga;| < ■ 

(3.15) 

Using Proposition 3.4 we can prove the following theorem establishing the 4>-system 
associated to the Lie algebra g. 

Theorem 3.6. Let q be a simple Lie algebra of ADE type, and let the solutions 
^(i) 

{x,E) : C —>■ i & I, have the asymptotic behavior (3.15). Then, the 

following identity holds: 

E) A'^^l\x, E)) = , for every i G L. (3.16) 

® 2 2 ' 
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Proof. Due to Proposition 3.4(b) and Theorem 2.4, the unique subdominant solu¬ 
tion to equation ( 2 . 1 ) in is 

(V-® +o(l)) , for X » 0 . 

Moreover, by equation (2.8), Corollary 2.7 and Proposition 3.4(c) we have that 

AU^%A^f+o{l)\ , for a; > 0 . 
2 2 V 2 2 / 

The proof follows by equation (3.14) and the uniqueness of the subdominant solu¬ 
tion. □ 

Example 3.7. In type An, equation (3.16) becomes 

^ 0 ^ for i = 1 ,..., n , 

^2 2 ' 

where we set = v]/("+i) = l. 

Example 3.8. In type D„, equation (3.16) becomes 

A^-f) for * = l,...,n-3, (3.17) 

mn-2 a 4'^””^^) = ® ^i("-l) 0 ^(") ^ (3.18) 

^2 2 ' 

A 4'^””^^) = ^'("-2) = TO„(^'^”| A ^-i”^) , (3.19) 

^2 2 ' ^ 22 ' 

where we set = 1 . 

Remark 3.9. The 4^-system in Example 3.8 is a complete set of relations among 

the subdominant solutions. The analogue system of equations obtained in [32, 

Eqs. (3.19,3.20,3.21)] lacks of the equation (3.18) for the node n — 2 of the Dynkin 

diagram. Note that equations [32, Eqs. (3.19,3.20)] are equivalent to (3.17) and 

(3.19) of the present paper. Indeed, it follows from the direct sum decomposition 3.1 

that rriiO L = ly(„- 2 ), for i = n — 1 , n, where t is the embedding ^ yM 

2 

(in fact, p("“ 2 ) ^ yr{n-i) ^ yrin) jjj particular case). The further identity [32, 
Eq. (3.21)[, which can be obtained from (A. 4), requires the introduction of another 
unknown and therefore equations [32, Eqs. (3.19,3.20,3.21)] are an incomplete 
system of n equations in n -f 1 unknowns. 

Example 3.10. In type E„, n = 6,7,8, equation (3.16) becomes 

A'i’f) = (g) , for i = l,...,n-4, 

^2 2 

^2 2 ' 

2 2 

A = ^-("- 3 ) 0 ^i(") , 

' 2 2 

^2 2 

where we set = 1 . 

4. The Q-system 

In this section, for any simple Lie algebra of ADE type, we define the generalized 
spectral determinants Q^'^\E-,I!.), i & I, and we prove that they satisfy the Bethe 
Ansatz equations, also known as Q-system. The spectral determinants are entire 
functions of the parameter E, and they are defined by the behavior of the functions 
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close to the Fuchsian singularity. More precisely, is the coeffi¬ 

cient of the most singular term of the expansion of (a:, E) in the neighborhood 
of cc = 0. 

In the case of the Lie algebra 5(2, the spectral determinant Q{E] t) was originally 
introduced in [14, 3], while the generalization to s[„ has been given in [10, 34]. An 
incomplete construction for Lie algebras of BCD type can be found in [9, 32]. The 
terminology generalized spectral determinants is motivated by the 5(2 case. Indeed, 
for the Lie algebra sb, the linear ODE (2.1) is equivalent to a Schrodinger equation 
with a polynomial potential and a centrifugal term, and the spectral determinant 
Q{E;i) vanishes at the eigenvalues of the Schrodinger operator. 

Before proving the main result of this section, we briefly review some well-known 
fundamental results from the theory of Fuchsian singularities. Here we follow ]22, 
Chap. III]. 

Remark 4.1. In this Section we assume that the potential p{x,E) = x^^'' is ana¬ 
lytic in a; = 0, namely Mh'^ G Z_|_. With this assumption, the point a; = 0 is simply 
a Fuchsian singularity of the equation. The case of a potential with a branch point 
at a; = 0 can formally be treated without any modification - see [3] for the case Ai 
- but the mathematical theory is less developed. 

4.1. Monodromy about the Fuchsian singularity. For any linear ODE with 
Fuchsian singularity at a; = 0, namely of the type 

+ H(x)) vlr(a;), (4.1) 

where A is a constant matrix and B{x) is a regular function, the monodromy 
operator is the endomorphism on the space of solutions of (4.1) which associates 
to any solution its analytic continuation around a small Jordan curve encircling 
X = 0. More concretely, if we fix a matrix solution Y{x) of the linear ODE (4.1) 
with non-vanishing determinant (i.e. a basis of solutions), then the monodromy 
matrix M is defined by the relation 

Y{e‘^^^x) = Y{x)M. 

Here, Y (e^’^'x) is the customary notation for the analytic continuation of Y (x) . 
We are interested in the invariant subspaces of the monodromy matrix M. By the 
general theory, the algebraic eigenvalues of the monodromy matrix M are of the 
form for any eigenvalue a of A. 

We said that an eigenvalue a of A is non-resonant if there exists no other eigen¬ 
value of a' such that a — a' G h. In this case the eigenvalue 6^^^“ has multiplicity 
one and there exists a solution of equation (4.1) which is an eigenvector of the 
monodromy matrix M with eigenvalue Such a solution has the form 

Xa(a:) =x“(xa + 0(x)) (4.2) 

for any Xa eigenvector of A with eigenvalue a. If the eigenvalues ai,..., are reso¬ 
nant, meaning that ai — aj G Z for j, j = 1..., fc, then the eigenvalue = ... = 

g27r2afc multiplicity k. The matrix M, in general, is not diagonalizable when 
restricted to the fc-dimensional invariant subspace corresponding to this eigenvalue. 
The corresponding solutions do not have the form (4.2) due to the appearance of 
logarithmic terms. 

We are interested in the linear ODE (2.1), which has a Fuchsian singularity with 
A = —i, £ G i). For any representation i G I, the eigenvalues for the action 
of £ on have the form A(£), where X G P is a. weight of the representation 
Note that for a generic choice of £ G\) and for distinct weights Ai and A2, the 
corresponding eigenvalues Xi{£), X 2 {£) are non-resonant. Indeed, generically, two 
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eigenvalues are resonant if and only if they coincide, and therefore they correspond 
to a weight of multiplicity bigger than one. 

4.2. The dominant term at cc = 0. Let w G , and let us denote by 

I* = max Re XU). 

Let also A* G P be such that I* = ReA*(£). If A* is a weight of multiplicity one 
in L{uj), then - by the discussion in Section 4.1 - the most singular behavior at 
X = 0 of a solution to equation (2.1) is x~^ . In this section, for any irreducible 
representation L{uj) of g and a generic £ € 1), we characterize the weight A* G P 
maximizing X{£) G Z. Moreover, we show that it has multiplicity one. These facts 
will be used to define the generalized spectral determinant Q^'^\E]£). 

It is well-known from the general theory of simple Lie algebras (see [19, Appendix 
D]), that given a generic element £ G 1), we can decompose the set of the roots R of 
g in two distinct and complementary sets (of positive and negative roots) 

R'l = {a G P I Re Q!(£) > 0} , R'^ = {a £ P | Rea{l) < 0} . 

Consequently, we can associate to this £ £ () a Weyl Chamber 211^, as well as an 
element W£ of the Weyl group (it is the element which maps the original Weyl cham¬ 
ber into We) and a (possibly new) set of simple roots = {we(ai) | i £ /}. The 
weight wi{uj) is the highest weight of L{ijj) with respect to the new set of simple 
roots A^, and due to the definition of P^, the action of any negative simple root 
wt{—oti), i £ I, decreases the value of Rewe{oj){£). Therefore, A* = we{uj) is the 
unique weight maximizing the function Re A(^), X £ P^j. This weight has multiplic¬ 
ity one since it lies in the Weyl orbit of w. 

In the case of a fundamental representation L[uji), i £ I, the weight we{u!i — at) 
belongs to P ^., and it has multiplicity one. Moreover, it is possible to show that all 
the weights in P^. different from we{u}i) are obtained from wi(uji — ai) by a repeated 
action of the negative simple roots we(—ai). We conclude that weiiOi — ae) is the 
unique weight maximizing ReA(^) on P^. \ {we{uii)}. We have thus proved the 
following result: 

Proposition 4.2. Let i £ 1) be a generic element, and let we the associated element 
of the Weyl group. Then the weight we(uj) £ Pu,, oj £ P'*", has multiplicity one and 

Rewe{X){£) >ReA(I'), 

for any weight X £ P^j, X we{u}). In the case of a fundamental weight to = oji, 
i £ I, the weight we{wi — ae) £ P^^ has multiplicity one and 

Rewe{uJi){£) > Rewe{uJi — ai){£) > ReA(I'), 
for any X £ P^j^, Xj£ wi{wi),wi{wi - at). 

Let G L{uji) be weight vectors corresponding to the the weights we{uJi) 

and we(uji — ai) respectively. As done in Lemmas 1.1 and 1.2 it is possible to show 
that Lijji) is a subrepresentation of /\ L{u!i) (respectively P(a;i)®^’^) and that 
A G L{r]i) (respectively £ Lijji)). Moreover, A and 

are highest weight vectors (with respect to A^) of the subrepresen¬ 
tation Lijji). Hence, we can identify these vectors (up to a constant) using the 
morphism mi defined in equation (1.4). From now on we fix the normalization of 
yb), ipU ), j G /, in such a way that 

A ■ (4-3) 
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(4.4) 


4.3. Decomposition of From the above discussion, and in particular form 

Proposition 4.2 and equation (4.2), we have that for any evaluation representa¬ 
tion i £ I and k £ there exist distinguished (and normalized) solutions 

X^k\x,E) and ip^^\x,E) of the linear ODE (2.1) such that they have the most 
singular behavior at a; = 0 and they are eigenvectors of the associated monodromy 
matrix. Explicitly, these solutions have the asymptotic expansion 

X«(a;,E;)=x-“‘(-^)W(x«+0(a:)), 

(^W(cc,e;) = + 0(x)). 

Since the parameter E appears linearly in the ODE (2.1) and the asymptotic be¬ 
havior at X = 0 does not depend on E, by the general theory [22] it follows that 
X^^\x, E;£) and ip^^\x,E-,(.) are entire functions of E. 

Using equation (1.17) and comparing the asymptotic behaviors we have that 

for every fc, fc' € C. For generic £ £ fj, the eigenvalues —wt{oJi){£), and —W£{uji — 
ai)(£) are non-resonant and therefore we can unambiguously define two functions 
Q^’'^E;£) and Q^'^\E]£) as follows. We write 

(x, E, £) = {E; £) x ^^^ (x, E) + (E; (x, E) + (x, E) , (4.5) 

where x*-®^(x, E) belongs to an invariant subspace of the monodromy matrix corre¬ 
sponding to lower weights. We call {E; £) and (E; £) the generalized spectral 
determinants of the equation (1.13). Since '^^'^\x,E,£), x^^'^{x,E) and (f^'''^{x,E) 
are entire functions of the parameter E, then also Q^'^\E;£) and Q^'^\E]£) are 
entire functions of the parameter E. 

Finally, the dt-system (3.16) implies some remarkable functional relations among 
the generalized spectral determinants. Indeed, we have the following result. 

Theorem 4.3. Let £ £ £) he a generic element. Then, the spectral determinants 
Q^'^{E-£) and QW(E;^) are entire functions of E and they satisfy the following 
QQ-system: 

Y[Q^^\e-,£)^'^ =oj-^^^Q^^\n-^E-£)Q^^\n^E-£) 

3^1 

-w^®-QW(D^F;;£)Q«(D-5F;;£), 

where 9i = we{ai){£ + h). 

Proof. We already showed that Q^'^\E;£) and Q^'^\E;£) are entire functions of the 
parameter E. By the definition of "i'^!\{x, E-,£) given in (2.8) and from equations 
(4.5) and (4.4) we have that 


(4.6) 


(x, E] £) = 


±^we{uji){e+h) f^{i} 


Q^^^{n^^-E-£)xf{x,E) 


(4.7) 


+ {p±hE- £)ipf (x, E) ■£ (x, E). 

2 

Here we used the fact that, since V^\ = V^\ then x^^{x,E) = x^\{x,E) and 

2 2 2 2 

ip^l\x,E) = {p^\(x, E). Note that v^'‘\x,E) belongs to an invariant subspace of 
2 2 

the monodromy matrix corresponding to lower weights. Equation (4.6) is then 
obtained by equating the component in in the i/j-system (3.16) by means of 
equations (4.5), (4.7) and the normalization (4.3). □ 
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We now derive the Q-system (Bethe Ansatz) associated to the Lie algebras of 
ADE type. Let us denote by G C any zero of Evaluating equation 

(4.6) at E = £1^2 E, we get the following relations 

1 ~ (4-8) 

Y[Q^^\n-2E;£f'^ =a;-^®^QW(L!-iE;,;^)QW(A,;£). 

jei 


Assuming that for generic £ &\) the functions (£’;£) and Q^'^\E;£) do not have 
common zeros, and taking the ratio of the two identities in (4.8), we get, for any 
zero Ei of Q^'^\E-,£), the Q-system 


JJ" fiPjCij 

i=i 


qU) (n^E*'j 


= -l, 


(4.9) 


where C = (C^)^ is the Cartan matrix of the Lie algebra g, and 




2Mh'^ ^ 

i£l 


{c-%e, = 


1 


2M/iV 


we{ujj}{£ + h ), 


j e / . 


Remark 4.4. The construction of the QQ system (4.6) works for any affine Kac- 
Moody algebra. However, formula (4.9) holds only in the ADE case. In order to 
define the analogue of equation (4.9) for the Cartan matrix of a non-simply laced 
algebra g, one has to consider, as suggested in [17], a connection with values in the 
Langlands dual of the affine Lie algebra g. 


4.4. The case £ = 0. The case of £ = 0 is not generic and thus the Q functions 
cannot be defined as in Section 4. It is however straightforward to generalize the 
generic case to ^ = 0, either directly or as the limit £ —>■ 0. Indeed, we can take 
the limit £ —^ 0 along sequences such that the element of the Weyl group wi = w 
is fixed. In this way, the limit solutions x^^\x,E,£ = E,£ = 0) of (2.1) 

satisfy the Cauchy problem 

(x, E,£ = 0)= vT , {x, E,£ = 0) = /f < . (4.10) 

Here u]" and /“' are respectively the highest weight vector of V^’'^ and the negative 
Chevalley generator corresponding to the element w of the Weyl group. Equations 
(4.10) define uniquely the two solutions. Then the spectral determinants 
are exactly the coefficients of {x = 0, E) with respect to the basis element 
v'^,f^Vi. Note that this implies a certain freedom in the choice of the spectral 
determinants In fact, the freedom in the choice of Q*^®^ corresponds to 

the elements of the orbit of uji under the Weyl group action, while the freedom 
in the choice of corresponds to the elements of the orbit of the ordered pair 
{uji,ai) under the same group action. 

For any choice of Wi = w, the functions Q^®\ satisfy (4.8) with £ = Q,we = w 
and the corresponding Bethe Ansatz equation (4.9), provided they do not have 
common zeros. 


5. g-AlRY FUNCTION 

We consider in more detail the special case of (2.1) with a linear potential 
p(x,E) = X and with £ = 0, which is particularly interesting because the sub¬ 
dominant solutions and the spectral determinants have an integral representation. 
Since in the 5 I 2 case the subdominant solution coincides with the Airy function, we 
call the solution obtained in the general case the g-Airy function. The classical Airy 
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function models the local behavior of the subdominant solution of the Schrodinger 
equation with a generic potential close to a turning point [16]; we expect the g-Airy 
function to play a similar role for equation (2.1). 


5.1. Q functions. In the special case of (2.1) with a linear potential p{x,E) = x 
and with .^ = 0, it is straightforward to compute the Q functions. By the discussion 
in Section 4.4, in order to define the Q functions we need to chose an element w of the 
Weyl group and thus the distinguished element vf of the basis of The spectral 
Q^'^[E) is then defined as the coefficient with respect to vf of 'I'(®i(a; = 0,E). Due 
to the special form of the potential we have 

4'(*i(x = 0,A) = 4'(*i(-£;,0). 


and thus, using formula (3.15), we get 


Q(*i(A) = '* C',(l + o(l)) as A^-oo, 


(5.1) 


where Ci is the coefficient of the vector \E\~^%j)^'''> with respect to the basis vector 
v'^. It is important to note that the asymptotic behavior of the functions is 
expressed via the eigenvalues which coincide with the masses of the (classical) 
affine Toda field theory, as we proved in Proposition 3.4. This is precisely the same 
behavior predicted for the vacuum eigenvalue Q of the corresponding Conformal 
Field Theory [31]. 

We conclude that the ODE/IM correspondence depends on the relation among 
the element A S g, the Perron-Frobenius eigenvalue of the incidence matrix B and 
the masses of the affine Toda field theory. 


Remark 5.1. Note that formula (5.1) does not coincide with the asymptotic behavior 
conjectured in [9, Eq 2.6] for a general potential. In fact, J. Suzuki communicated 
us [33] that the latter formula is conjectured to hold only for potentials such that 
< 1, in which case the Hadamard’s product [9, Eq 2.7] converges. Clearly in 
case of a linear potential ^ > 1. 

Remark 5.2. A connection between generalized Airy equations and integrable sys¬ 
tems appears, among other places, in [23], where the the orbit of the KdV topolog¬ 
ical tau-function in the Sato Grassmannian is described. We note that an operator 
of the form (1.13) - specialized to the case of algebras of type and for I = —h/h'^, 
M = 1/K^ - also appears in that paper. It would be interesting to obtain a more 
precise relation between these results and those of the present paper. 

5.2. Integral Representation. Given an evaluation representation of g, we look 
for solutions of the equation 

4''(x)-I-(e-I-xeo)'I'(a:) = 0 , (5.2) 

in the form 

^{x) = f e-^^^is)ds, (5.3) 

J C 

where $(s) is an analytic function and c is some path in the complex plane. Dif¬ 
ferentiating and integrating by parts we get 



where the last term is the evaluation of the integrand at the end points of the path. 
We are therefore led to the study of the simpler linear equation 

( - s-I-e-I-eo^)$(s) = 0, (5.5) 

which we solve below for two important examples. 
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5.3. An-i in the representation In this case equation (5.2) is already 

known in the literature as the n-Airy equation [16]. Using the explicit form of 
Chevalley generators Ci, i = 0,... ,n — 1, which is provided in Appendix A.l, and 
denoting by $i(s) the component of ^(s) in the standard basis of C", we find the 
general solution of (5.5) as 

s^+l 

$i(s) = Ke"+i, $j(s) = $i(s), i = 2,...,n, 


where k is an arbitrary complex number. If we choose the path c as the one 
connecting e~"+i oo with e"+i oo, then the integral formula (5.3) is well-behaved, 
and the boundary terms in (5.4) vanish. Thus, the s[„-Airy function has the integral 
representation 


n + l oo 




27ri rT+i 




n + 1 

-*+1 ds. 


j = l,...,n. 


(5.6) 


+ 1 OO 


In case n = 2, the latter definition coincides with the definition of the standard 
Airy function. By the method of steepest descent we get the following asymptotic 
expansion of T for x ^ 0: 





j = l,...,n. 


(5.7) 


Due to Theorem 2.4, the s[„-Airy function coincides with the fundamental solution 
of the linear ODE (5.2). In fact the asymptotic behavior of the Airy function 
coincides with the asymptotic behavior of the function (for this computation 
one needs to use the explicit formula for h which is given in equation (A.l) ). The 
solutions k € h, are obtained by integrating (5.3) along the contour obtained 

2k-Ki 

rotating c by e ’*+i. 


5.4. Dn in the standard representation. The second example is the D„-Airy 
function in the representation Using the Chevalley generators e^, i = 0,..., n, 
as described in Appendix A.2, and denoting the jth-component of $ in the 
standard basis of C^", we find that the general solution of (5.5) reads 

s2n-l 

$1 = K 2 e 2n-l 

= s-’“^$i(s), j < n - 1, 

n-l 

= -^4*1(3), $„+! = S"“l 4 >i(s), 

= s"+-^“^4>i(s), j<n-l, 

/ 2«-2 I \ 

(s). 

where At S C is an arbitrary constant. Let c be the path connecting e~ 00 with 
00. Then along c the integral formula (5.3) is well-behaved and it thus defines 
a solution of (5.2) that we call ZI„-Airy function. Moreover, the boundary terms 
in (5.4) vanish. As in the case A„, the standard method of steepest descent shows 
that such solution is exactly the fundamental solution to equation (5.2). The 
solutions fc S Z, are obtained by integrating (5.3) along the contour obtained 
rotating c by . 


Appendix A. Action of A in the fundamental representations 

In this appendix we give an explicit description of the maximal eigenvalues of A 
in the fundamental representations and of the corresponding eigenvector. 
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A.l. The An case. The simple Lie algebra of type T„, n > 1, can be realized as 
the algebra of (n + 1) x (n + 1) traceless matrices 

0 = slra+i = {A G Mat„+i(C) I Tr A = 0} , 

where the Lie bracket is the usual commutator of matrices. The dual Coxeter 
number of g is /i^ = n + 1. Let us consider the following Chevalley generators of g 
{iGl = . ,n}): 

fi — 1 hi — Eli £'2+1,Z+1 ; — -^2,2 + 1 ; 

where Eij denotes the elementary matrix with 1 in position (i, j) and 0 elsewhere. 
It is well-known that the representation L{uji) is given by the natural action of g 
on L(a;i) = Moreover, we have that 

i 

L{iOi) = I\L{loi) , i Gl. 

We denote by Uj, j = 1,... , n-|- 1, the standard basis of and hy Vi, i G I^ the 

highest weight vector of the representation L(uji). Then, we have: 


Vi = ui A U2 A ■■■ A Ui. 

The set of Chevalley generators for g is obtained by adding to the Chevalley gen¬ 
erators of g the following elements: 

/o = Ei^n+lt ^ , ho = 2c— Ell + En+l^n+l , Cq = En+l^lt . 

The element h G i) satisfying relations (1.14) is 


u A- I 

h = drag ( - -,- 


n — 1 n, 

2 ’ V ■ 


(A.l) 


Recall that A = eo+ei + - ■ ■+en and, from Example 1.3, that = /\* L(a;i) , 

~T~ 

for i G L In particular, and we set 

n+1 

^(1) =J2ujG . (A.2) 




Then, it is easy to check that . By Proposition 3.4(a), G is 

the unique (up to a constant) eigenvector corresponding to the maximal eigenvalue 
= 1. Furthermore, using equation (3.10), we can check that for every i G I, 


A« = 


sm 


(^) 


sm 


(n+l) 


(A.3) 


is a maximal eigenvalue of A in the representation Using equations (3.3) and 
(3.2), the corresponding eigenvector is easily checked to be 

A A • • • A A G . 


A.2. The Dn case. Let n G and consider the involution on the set {1,..., 2n} 
defined hy i i' = 2n + 1 — i. Given a matrix A = (Ay)^"_^ G Mat2ri(C) we 
define its anti-transpose (the transpose with respect to the antidiagonal) by 


^at 


(^) 


2n 

i,j=l ’ 


where A)®* = Ajfi' . 


S = {Ekk+Ek'k') ■ 

k^l 
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Let US set 





Following [15], the simple Lie algebra of type can be realized as the algebra 

0 = 02„ = {A e Mat2„(C) I AS + = 0 } , 

where the Lie bracket is the usual commutator of matrices. The dual Coxeter 
number of g is = 2n — 2. For i,j € I = {1,2,, n}, we define 

+ {-ly+^+^E,,,, , F,, = , 

and we consider the following Chevalley generators of 0 : 

fi — Eij,.\i , hi — Fa 2+1 ? 2+1 j f 1 , 

fn — 2^22/ (21 — 1)^ ) l,n —1 “t“ Eji ji , €ji — . 

It is well-known that the representation L{uji) is given by the natural action of 0 
on L{uji) = Moreover we have that 

i 

L{ujy) = l\L{uJi), i = l,...,n-2, 

while L{ojn-i) and Liujn) are the so-called half-spin representations of 0 . We denote 
by Mj, j = 1 , • ■ •, 2 n, the standard basis of and hy Vi, i € I , the highest weight 
vector of the representation L(uji). It is well-known that 


Vi = Ui A U 2 A ■■■ A Ui, i = 1,... ,n — 2 . 

The set of Chevalley generators for 0 is obtained by adding to the Chevalley 
generators of 0 the following elements: 

/o = 2 F( 2 ,i)'_( 22 i-l)'t ^5 ho = 2c— Ell — F 22 , Cq = —F2n-l,2nt ■ 

Recall that A = Cq -I- ei -I- • • • -I- Cn and, from Example 1.3, that Rb) = /y* , 

~T~ 

for i = 1,... ,n — 2, = L(ujn-i)^. and = Liyjn)^ are the evaluation 

representations at C = (—1)” of the half-spin representations. In particular, = 
C^”, and we set 

n — 1 ^ 

^ {Uj + Un+j) + 2 (■“" + “2n) ■ 

i=i 

Then, it is easy to check that By Proposition 3.4(a), is 

the unique (up to a constant) eigenvector corresponding to the maximal eigenvalue 
= 1. Furthermore, using equation (3.10), we can check that for every i = 
1 ,... ,n - 2 , 



is the maximal eigenvalue of A in the representation and the corresponding 
eigenvector is easily checked to be 

•0b) = 0b)_^ /\ 0b]_^ A • • • A 0^ A 0^ G pb) ^ 

2 2 2 2 

as follows from equations (3.3) and (3.2). For the fundamental representations 
P("- 1 ) and P(") (see [19] for the definition), it is easy to check, using again equation 
(3.10), that 

\("-i) _ \(n) __ ^ 

" ^2si„(++)' 

Finally, we note that there is another important irreducible representation for the 
Lie algebra 0 , which we denote by p 1^^ It is not hard to show 
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that is the maximal eigenvalue for the action of A on Its 

corresponding eigenvector is 

^(n-l) _ A A • • • A A • 

2 2 2 2 

Moreover, the following isomorphism of representation (see [19]) 

(g) T/(”) ^ [/("-!) 0 y("-3) 0 ... ^ 

implies that there exists a morphism of representations m such that 

TO ■ (A.4) 

A.3. The E case. In Table 2 we give the explicit form of the eigenvalues A*^*^’s, 
obtained from (3.10). The corresponding eigenvectors are computed in a GAP [21] 
file available upon request to the authors. 



Ee 

Er 

Es 

A(i) 

1 

1 

1 

A(2) 

sin(f) 

'^“(fi) 

sin(f ) 


A(3) 

sin(f ) 

'^“(fi) 

sin(f) 

sin(,!^) 

sin(^) 

A(4) 

sin(f ) 
sin(^) 

sin(^) 

sin(^) 

"“(if) 

"“(if) 

A(5) 

sin(f) 

sin(^) 

sin(^) 
sin(f ) 

"“(f) 

"“(is) 

A(6) 

1 

sin(f)sin(^) 

sin(f)sin(f^) 

"“(f) 

"“(ff) 

A(7) 


sin(%^) 

sin(f) 

sin(f )sin(f^) 
sin(f^)sin(fj) 

A(8) 



"“(f) 

"“(fs) 


Table 2: Maximal eigenvalues for simple Lie algebras of E type 
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